Macroscopic view of light pressure on continuous medium 
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The ambiguity of macroscopic description of light pressure on continuous medium originates from 
the uncertainty of dividing the energy-momentum tensor of electromagnetically excited matter into 
material and field parts or, equivalently, the total acting force into pressure and deformation terms. 
We show that although there exists a continuum of formally correct formulations, one can adopt the 
appropriate form of the macroscopic field stress tensor that allows unified description of pressure 
during elementary light-matter interactions, such as reflection/refraction, absorption and nonlinear 
conversion. The proposed expressions for the pressure force are simple, convenient and compatible 
with the polariton momentum ftk. The corresponding electromagnetic momentum density (|14[1 
generalizes Rytov's definition for right-handed and left-handed frequency dispersive media. 

PACS numbers: 42.50.Wk; 42.25.Bs; 78.20.Bh 



Since the pioneering works by Lebedev the light 
pressure has been the most substantial observable con- 
sequence of the momentum transfer by electromagnetic 
(EM) waves in vacuum. For decades afterwards, the phe- 
nomenon was supposed to be primarily of general im- 
portance with uncertain practical value reduced to e.g. 
gigantic solar-sail propelled spaceships [2j. Later it be- 
came clear that observable and useful manifestations of 
the light pressure occur actually on the submillimeter 
scale, where powerful laser beams can be used as tweez- 
ers for microparticles Nowadays, the optomechanics 
has developed into a broad area with numerous potential 
applications in nanotechnology (see e.g. Q). 

In the same way as for 'pure' photons in vacuum, the 
concept of light pressure in continuous medium has to 
rely on momentum carried by the corresponding quasi- 
particles - polaritons. However, the definition of the po- 
lariton momentum remains the subject of ongoing dis- 
cussions. The crucial obstacle is the lack of established 
unique definitions of momentum density and flux of EM- 
field in continuous medium. 

Over the years, various definitions have been proposed 
and analyzed [1-15[. Two simplest and most prominent 
variants belong to Minkovski and Abraham. Accord- 
ing to Abraham [eJ , the momentum density of EM- wave 
reads Ga = E x H/(47tcV, while Minkovski's alternative 
is Gm = D x B/(47rc) Q. The latter has been gener- 
alized by Rytov [7] for substantial frequency dispersion 
involving terms with frequency derivatives of permeabil- 
ity and permittivity (see also More recently, other 
less trivial variants have been proposed (e.g., a half-sum 
of Ga and Gm 0) although it has been argued that 
only Ga and G m possess special physical meanings and 
other 'rival forms' are excessive [III ]. 

Surprisingly, no direct controversy seems to follow from 
any reasonable definition, which supports the common 
belief that "the preferred form is therefore effectively a 
matter of personal choice"\l'^. To avoid this uncertainty 
when solving practicalproblems (e.g. calculating the 
forces in confined fields [16j ) one often encloses a finite ob- 



ject (microparticle) by a surface in vacuum and uses the 
vacuum relations to obtain the overall force and torque. 
Then, however, the details of the momentum transfer re- 
main unrevealed and one cannot trace the contributions 
arising from reflection/refraction at the surfaces, absorp- 
tion in the volume or finite pulse duration or explain the 
roles of key parameters (particle shape, refractive index, 
absorption, etc.). 

In some applied domains the above unclarity has not 
even been fully recognized. Thus throughout decades 
in the literature on nonlinear optics the phase matching 
condition is mixed with the light momentum conserva- 



tion law (see e.g. the classical (17J, |18|, or the recent 
(ljj). which imposes the polariton momentum directly 
expressed by the polariton wave vector as Kk. However, 
the link between phase matching and momentum transfer 
during nonlinear processes has not been established. 

The topic has got recently another intriguing dimen- 
sion with emergence of metamaterials. While in conven- 
tional dielectrics the ambiguity of polariton momentum 
can be seen as a minor quantitative issue, in the left- 
handed media (LHM) with negative permeability and 
permittivity, the vectors of polariton phase and group 
velocities are antiparallel 0, [2(| and it becomes a matter 
of momentum and pressure signs || [2l[ |22| ■ 

In this Letter our aims are: (i) to analyze the actual 
degree of freedom in defining light momentum and pres- 
sure in continuous medium and to show that there exists 
a continuum of possibilities of formally self-consistent de- 
scription fulfilling the necessary conservation laws; (ii) to 
propose a unified approach that yields simple and conve- 
nient expressions for the light pressure during elementary 
light-matter interactions. In this approach, the forces 
acting on inhomogeneous medium, as well as during light 
reflection/refraction and absorption, pulse propagation 
and nonlinear conversion can all be presented by simple 
self-consistent expressions that are compatible with the 
polariton momentum Kk and the EM-field momentum 
density as simple as Eq. ([f"4|) applicable to both right- 



handed and left-handed frequency dispersive media. 
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Several relations on the microscopic level form the ba- 
sis for the following macroscopic analysis. Thus the mi- 
croscopic forces acting on matter can be conveniently de- 
rived from the Maxwell's stress tensor 12311 : 



Uj = (An)- 1 [e.ej + bibj - (e 2 + b 2 ) /2] 



(1) 



expressed by the microscopic electric and magnetic fields 
e(r, t) and b(r, t). The stress tensor equates with the mo- 
mentum flux tensor taken with the opposite sign, which 
can be illustrated by considering (JXJ) for a monochro- 
matic plane wave in vacuum (wave vector k — lj/c, fields 
e = — c/oj k x b and b = c/lu k x e): 



,pw 



(47r) *(e x b)ikj/k = -g%Vj 



(2) 



Here v = ck/fc is the group velocity and g=ex b/(47rc) 
is the microscopic density of momentum. The latter form 
presents the stress tensor as a dyadic product of the mo- 
mentum density with the velocity of its transmission, i.e. 
as the momentum flux tensor. 

To define the momentum quant in this otherwise classi- 
cal picture one introduces the energy quant Huj and writes 
the energy flux and the momentum flux as proportional 
to the photon number flux N: S = e x b c/4ir — Nhiu, 
and tij = —NiPj. Comparing with Eq. @ yields the 
quant momentum p = ftk. 

The divergence of Eq.([T]) yields the microscopic mo- 
mentum balance 12311 : 
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(3) 



where i"l is the microscopic Lorentz force density acting 
on microscopic charge and current densities p nllc and j m i C : 



h = Pmi C e + (jmic x b)/c. 



(4) 



Similarly, to derive the microscopic angular momen- 
tum density and torque one introduces the angular mo- 



mentum tensor 



IE Ha Hi: 



(5) 



where 5ijk is the antisymmetric Levi-Civita symbol. The 
microscopic angular momentum balance reads 
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(6) 



where L = r x i"l is the microscopic torque and 1 = r x g 
is the angular momentum density. 

For the EM-fields periodic in time, the last terms in 
Eqs. ([3j [6]) vanish upon time averaging and the average 
total force and torque acting on a finite body in vacuum 
can be evaluated from the stress tensor distribution out- 
side the body as: 



F|,„ = I dV t L = i t-</S 



(7) 



= / dVcj> L = 
Jv 



m • dS, 



(8) 



respectively, with the integrals taken over the volume V 
including the body and the surface S in vacuum enclosing 
the body, while the bar stands for the time averaging. 

Macroscopic theory of light pressure has to comprise 
macroscopic analogues of Eqs. © and ([5]) with appro- 
priate macroscopic force density. Unfortunately, a gen- 
eral macroscopic averaging of the above is impossible 
due to the quadratic dependence of the quantities on 
microscopic fields and densities of charge and current. 
Note that such an averaging can be performed within 
particular microscopic models of the medium (see e.g. 
|l4 l25l . [26 1) but yields different model-specific results. 
This is not surprising since the procedure accounts for 
all forces arising in a body subjected to EM-fields: those 
responsible for the momentum exchange and those that 
tend to deform the body without momenta transfer. 

The possibility of various formally self-consistent ap- 
proaches is recognized to be the consequence of the am- 
biguity of separating the energy-momentum tensors of 
EM-field and matter 12, Uj]. It is only the total tensor 
of the whole system that has a clear physical meaning. 
The separation line is uncertain but we find that only 
few general restrictions are sufficient to ensure the fulfill- 
ment of conservation laws. Namely, the coordinate part 
of the energy-momentum tensor - the EM stress tensor 
(macroscopic analogue of Maxwell's tensor (J])) - must: 
(i) be quadratic in the macroscopic fields E, B, D, H; (ii) 
tend to Eq. ([TJ in vacuum; (iii) be symmetric. 

For any such a bilinear form T(E,B,D,H) one can 
write the macroscopic momentum balance: 



(9) 



which serves then as definition of the pressure force Fl- 
For a finite body enclosed by a surface S, T = t on S 
(in vacuum), and the total force calculated as the in- 
tegral J v F^dV will always be equal to F to t given by 
Eq.([7]) thus fulfilling the total momentum conservation 
automatically. 

As usual in the field theory [24( , the conservation of the 
total angular momentum will also hold true provided that 
is symmetric. Indeed, for the macroscopic angular 
momentum tensor My = 5ik n ZkT n j, one can introduce 
the macroscopic torque $l via: 

d 



dx 



Mi. 



(10) 



which equals the torque exerted by the pressure force, 
$ = rxFi, and sums up to the total torque J v &lcIV 
that automatically equals Eq. ((5J since M = m on the 
surface S in vacuum. 

Using the macroscopic force and torque defined from 
Eqs. © and (fTU|) one can consider elementary light- 
matter interactions, introduce an appropriate EM-wave 
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momentum density to ensure the momentum conserva- 
tion, etc. Obviously, there can be as many different for- 
mulations as there are possibilities to write the appropri- 
ate bilinear form T(E, B,D,H), i.e., an infinite contin- 
uum. Existing theories as those in Refs. d,Gil 15.l27l.l28j 
can all be obtained using this strategy. In this sense, they 
all are equally correct: however different spatial distri- 
butions of light pressure and torque they yield, the total 
force and torque acting on the whole body is the same. 

For example, choosing the stress tensor 
[EiEj + BiBj - 6ij/2 (E 2 + B 2 )] /4tt yields the 
macroscopic pressure force pE + (j x B)/c, as if the 
macroscopic fields act on the macroscopic densities of 
bound charges p and currents j. Essentially, it is the 
similarity with (HJ) that has motivated the authors of 
Refs. jlfl 121 l28( to use such approach. However, it 
cannot really be obtained by macroscopic averaging 
of ((4]) and its severe drawbacks become obvious as 
one applies it to particular light-matter interactions. 
Thus the force is nonzero already for a superposition of 
EM-waves of the same frequency propagating through 
lossless homogeneous medium in different directions, 
i.e., the momentum exchange occurs in the absence of 
real interaction. Such 'virtual' terms must be then taken 
into account to ensure the momentum conservation and 
calculations become cumbersome. 

A much clearer picture follows from the stress tensor : 



BiBj- 



Tij = (87T)- 1 [EiDj + DiEj + B Z B 3 

^■(E-D + H-B)], (11) 

first proposed by Rytov Q- Similar but unsy mmetrized 
versions have been considered e.g. in Refs.ll2l andl29l and 
for monochromatic fields in homogeneous media (|11[) re- 



duces to the classical form [13J, [16j, |23j , but it is the ex- 
pression (fTTj) that yields correctly all the relations below. 

Consider first Eq. (|11[) for a monochromatic transver- 
sal plane wave traveling through an isotropic medium 
with e and /i of the same sign, i.e., with the wavevec- 
tor k — sfeji ut/c and the fields D = —c/lj k x H, 
B = c/lo k x E, D = eE, B = /xH. The stress tensor 
(fTTj) then reads 

TFT = ~—ME x H)j = -piNj, (12) 



4irui 



where the polariton density flux N = S/(hu) is expressed 
by the Poynting vector S, and the latter form is the diadic 
representation with the polariton momentum p = fik. 

On the other hand, the polariton dispersion law k 2 = 
£fi to 2 j c 2 determines the group velocity 

-l 



duj(k) 
3k 



2c 2 



d(uje) d(u)fj.) 



(13) 



and one can write (H 



-GiVj (compare with 



Eq. ((2J) with the momentum density of the EM wave 



G = (E x H) 

8vrc ' 



d(tue) d{u>n) 



du> 



du> 



(14) 



This most general form reduces to Rytov's definition Q 
in ordinary media (for positive e and /i) and further to 
Minkovski's G M upon neglecting the dispersion. 

Next, taking space derivatives of (ITT1) it is easy to 
obtain the macroscopic differential force-momentum bal- 
ance with the total force 

F = (8tt)" 1 [D(V • E) + B(V ■ H) - E x curlD 

D x curlE H x curlB B x curlH]. (15) 

For an arbitrary combination of monochromatic 
fields in a homogeneous isotropic transparent medium, 
this force reduces to the full time derivative F = 
(Anc)" 1 d/dt(D x B) and its time -average F vanishes 
identically. Therefore, no virtual forces arise for a super- 
position of non-interacting EM-waves. 

In a transparent and isotropic but inhomogeneous 
medium, 



1 (e 2 \7s + B 2 Vfi 



(16) 



i.e., the averaged force always points opposite the di- 
rection of the gradients of fi and e, and, in particular, 
normally to interfaces. 

To obtain the force acting on an abrupt xy-plane in- 
terface between two media, one can introduce a transient 
layer of infinitesimal thickness S with e(z) and fi(z) con- 
tinuously varying from £\ and [i\ at z < —5/2 to £2 and 
fi2 at z > 5/2 respectively. Using the tangential field 
components E r and H T and normal induction compo- 
nents D z and B z that stay constant throughout a thin 
layer, the normal force can be expressed as 

F int = j F z dz= — [(ei - e 2 )£? + (Mi - m)H 2 + 



(5/2 



8tt 



(l/e 2 - l/ 6l ) D 2 + (1/ M2 - 1/ M1 ) B 2 }. (17) 

It does not depend on the particular structure of the 
transient layer and can be evaluated, for example, for 
the fields arising when an arbitrarily polarized EM-wave 
is incident on the interface at an arbitrary angle. Sim- 
ple calculations show that this force coincides with that 
evaluated from the simple 'corpuscular' relation 



Kk N 0z - Kk t N tz - Kk r N r 



(18) 



where ko, r ,t are the wavevectors of the incident, reflected 
and transmitted waves respectively, and ATo,r,t z are the 
^-components of the corresponding polariton number 
fluxes. 

In a weakly lossy homogeneous medium, i.e., for fi = 
fi' + ip" and e = s' + is" with small imaginary parts, 
the complex representations of the oscillating fields: E = 
(E e- itJt + c.c.)/2, D = (eE e- lwt + c.c.)/2, etc., substi- 
tuted in Eq. flTBl yield after time-averaging 



4-7TC 



Ve")(E x H). 



(19) 
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This force is always collinear with the average momentum 
density (H~4|) and corresponds to the absorbed polariton 
momentum. Thus for a linearly polarized monochromatic 
plane wave it is proportional to the absorbed energy Q = 
uj (47r) ~ 1 (e"~EP + as F = Qk/w, which means that 

the absorption of energy quant Huj is accompanied by the 
absorption of momentum quant fik. 

For a finite slowly modulated pulse of monochromatic 
light propagating through isotropic non-lossy medium, 
one introduces complex slow varying amplitudes, e.g., the 
electric field E = l/2[E (i)e -iu * + c.c.]. The frequency 
dispersion modifies the relations between field amplitudes 
and inductions [IH so that the slow-varying electric in- 
duction amplitude reads D (t) = eE (£) +i de/duj E (t) 
and the slow amplitude of its first time derivative equals 
(dD/dt) (t) = -iweE (t) + d(oje)/duj E (f). Substitut- 
ing these amplitudes together with the similar magnetic 
ones into Eq.Q yields a very simple result: 



dG 

at ' 



(20) 



where the time averaged momentum density G is given 
by Eq. (|14p taken with the slow varying field amplitudes. 
This represents a generalization of the Abraham force. 

Similarly, one can evaluate the momentum transfer 
during light-light interactions in nonlinear media. For 
simplicity, consider interaction of collinear waves prop- 
agating along z-axis and polarized linearly along x-axis 
and assume that the energy exchange is dominated by 
a single channel for which the set of frequencies {w„} 
fulfills the condition ^2 n l n uj n — with integer l n and 
the wavevectors are weakly mismatched: l n k n = q. 
For definiteness, consider the dielectric nonlinearity with 
nonlinear dielectric susceptibility ^-( A/ ) of the order M. 

Then the electric field can be presented via the 
slow-varying with z amplitudes: E = ^ n E n = 
l/2J2 n [ E no(zy knZ ~ iuJnt + cc] which obey the system 
of equations: 



dE n0 27Ti0J l„ (l , AP NL 



(21) 



Here the slow-varying amplitude of nonlin- 
ear polarization at frequency Lu n reads f^o = 
2 -M e - Bg n( yi? , x (M) (£f(o) i-|y nn; ^ (£n/o) - ssn (y in , ) 

the order of the nonlinearity is M = J2 n I'nl — 1> f° r ne §" 
ative field powers the complex conjugate amplitudes E* 
must be substituted, and Xn^ is the ir-component of the 
susceptibility tensor taken at appropriate frequencies: 

AM) _ (M) (.,.{., ,\\ 
Xn — Xx... x\Un, \ u n'JJ- 

Upon substituting D = J2 n e(w„)£ n +47rP NL and B = 
J2 n n{<jJ n )Hn into the force (|T5|) only the terms 



AP 



NL 



2 V dz 



dz 



contribute to the time average. Neglecting higher pow- 
ers of nonlinearity one can conveniently express the force 
density by the polariton number density fluxes, N n — 
S n /(hu) = c/(8irhu n )y/e(uj n )/n(uj n )\E n0 \ 2 : 



i? Y^fc/ dNn n \ dN > 



dz 



(23) 



(22) 



For the perfect phase matching, q — 0, only the first 
term is nonzero and the force has intuitively clear 'cor- 
puscular' form compensating the possible momentum un- 
balance during conversion and vanishing, for example, for 
phase-matched M-th harmonic generation. In presence 
of phase mismatch, however, the second term in Eq. (|2"3")l 
gives a nonzero contribution as the medium acquires the 
mismatched momentum. This establishes the link be- 
tween the phase matching and the momentum conserva- 
tion, which has been intuitively adopted e.g. in Refs. 
IT9I and numerous others. 

Finally, consider the implications of the above for 
LHM. First, the momentum density (fl"4|) is parallel to 
the Pointing vector in ordinary media and antiparallel in 
LHM (e < 0, fj, < 0). Indeed, the involved frequency 
derivatives must stay positive to ensure the stability of 
the EM energy |18J, and, therefore, the sign of the square 
brackets follows that of e and /1. Next, the force acting 
on an interface between media of different handedness ac- 
cording to (|T7|) always points in the direction of the LHM. 
For a slab of lossless LHM, however, the forces from the 
input and output interfaces partially compensate to en- 
sure the total momentum conservation. In a lossy LHM 
the force (|T9)) is opposite to the Poynting vector. For 
phase matched harmonic generation in non-dissipative 
LHM, no pressure is exerted in the bulk although an 
actual reversal of the energy flux can take place if the 
medium is left-handed not for all waves. In the latter 
case for the material acting as a nonlinear mirror (3pj 
the total light pressure arises solely at the surfaces in 
accordance with Eq. (fT8|) at each frequency. 

In conclusion, we have shown that although there ex- 
ists a continuum of possibilities to describe the momen- 
tum exchange between EM radiation and continuous 
medium, one can formulate an approach that provides 
a simple and intuitively clear picture. It is compatible 
with the polariton momentum Kk, yields the compact 
general form (|14j) of the momentum density of EM wave 
and allows simple representation of light pressure dur- 
ing reflection/refraction, absorption and nonlinear con- 
version. This enables one, in particular, to divide the to- 
tal momentum transfer in a complex real situation into 
elementary processes yielding detailed understanding of 
the roles of key parameters and conditions. 

We are grateful to our colleagues V.E. Dmitricnko, B.I. 
Sturman, and V.A. Chizhikov for stimulating discussions, 
criticism and helpful advices. 
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